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SUMMARY 


Two programs for the cumulative binomial function 


n 


S —si— P 1 
Z— / j!(n - j)! 


p 3 (l - p) n “ 3 


]=r 


are presented. The direct program computes C for specified values of n, r, and p. 
The indirect program computes p for specified values of n, r, and C. The programs, 
which are suitable for n :£ 1000, are useful in applications such as reliability, industrial 
sampling, and statistics. 


INTRODUCTION 

Calculation of probabilities using the cumulative binomial function is useful in such 
applications as reliability, industrial sampling, statistical-hypothesis testing, and 
confidence- interval determination. The cumulative binomial function is 


C = 



p 3 (l - p) n " 3 


Two calculator programs are presented for the occasional user of the binomial func- 
tion who has available a Hewlett-Packard 9100 A desk-top calculator. The direct pro- 



gram computes the cumulative binomial sum C for given values of n, r, and p. The 
indirect program uses an iterative procedure to determine the value of the binomial pa- 
rameter p for specified values of n, r, and C. 

Instructions for use of the programs are given along with a discussion of limitations 
and applications. A more thorough treatment of the indirect problem is presented in the 
reference along with extensive tables of p and several applications. 


SYMBOLS 

C cumulative binomial probability 

cumulative binomial sum using p^ as the binomial parameter 
j summation variable 

n upper limit on binomial summation 

p binomial parameter; element success probability 

r lower limit on binomial summation 

Subscripts: 
f final value 

i iteration number 

0 initial value 


THEORY 


The probability C of at least r successes among n independent identical ele- 
ments of a system is given by the cumulative binomial function 


where 



0 


( 1 ) 


( 2 ) 


2 



and p is the element success probability. The negative binomial distribution is equiva- 
lent to equation (1); its equation is 



Equation (3) is used in both the direct and indirect programs because it requires fewer 
program steps than equation (1). In the expansion of the summation, the first term is 
always unity. Successive terms of the surfimation are obtained by multiplication by 

(1 - p) (4) 

i 

DIRECT PROGRAM 

The function of the direct program is to determine the value of C for given values of 
n, r, and p by summing the terms of equation (3). The user instructions and program 
steps are described in figure 1. The three inputs n, r, and p are loaded by the user 
into the z, y, and x registers, respectively, before beginning execution. After load- 
ing, the user must press "Go to (0)(0)" and "Continue. " The display is C in the x 
register and zeros in y and z. 

The following sample cases illustrate the use of the direct program. The times 
given are the approximate delays between pressing "Continue" and the display of the 
result. 


Case 1: n = 10 

C = 0.967 206 502 

r = 6 


p = 0.8 

less than 1 second 

Case 2: n = 100 

C = 0.999 192 426 

r = 80 


p = 0. 9 

about 2 seconds 

Case 3: n = 1000 

C = 0. 526 599 080 

r = 900 


p = 0. 9 

about 6 seconds 
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To 4-0 
From 3 - 


Enter program: (Starting address is 0 - 0) 
Enter data: n-^z, r*^y, p-^-x 
Press: Go to (OHO) or End 
Press: Continue 
Display: 

0 z 

0 y 

C x 

(a) User instructions. 
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(b) Program steps. 

Figure 1. - Direct-program user instructions and program steps. 










Limitations 


The direct program incorrectly sets C to zero whenever p r is less than 10"^, 
the limit of the calculator. This occurs because C is proportional to p r (see eq. (3)). 
Figure 2 shows the allowable pairs of p and r for r from 10 to 1000. The curve 
separating the allowable and unallowable regions is the equation p r = 10"®®. When the 
calculator displays either C = 0 or C = 1, the input data should be checked. 



Figure 2. - Maximum r and minimum p values allowed in direct program. A ilowable values determined from 
p r >lQ"". 


Error Light 

The error light comes on only when the input data are improper. This occurs when 
r is greater than n or when p is outside the allowable range of 0 to 1. 


INDIRECT PROGRAM 

The user instructions and steps of the indirect program are presented in figure 3. 
This program determines the value of p for given values of n, r, and C. This is done 
with the aid of the Newton-Raphson iterative procedure described in the reference. The 

5 




Enter program*. (Starting address is 0 - 0) 
Enter data: n-^z, C-»-x 

Press: Go to (OKO) or End 
Press: Continue 

Pause display: 

(C-Cf) *• 2 

ipi+i-pi 1 — y 

p, X 

Final display: 

(C - 6 h > z 

I Pf ~ Pf-1 1 y 

P f x 

(a) User instructions. 



(b) Program steps. 

Figure 3. - Indirect-program user instructions and program steps. 
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(b) Concluded. 
Figure 3. - Concluded. 



















(i + I)** 1 estimate of p is determined from the i^ 1 estimate of p as follows: 


C - C 4 


P i+ l = Pi + 


aPi 


( 5 ) 


where 


and 


■■•■I 

i=o 


( 6 ) 


5"CK' 


(7) 


In order for convergence to be assured, the initial guess for p is chosen at the in- 
flection point on a curve of C as a function of p; that is, 


?0 


r - 1 
n - 1 


( 8 ) 


Equation (5) is then used to determine all subsequent approximations of p until the in- 
crement of p is less than 10" ; that is, 


|P i+1 -Pil <10' 6 (9) 

The three inputs n, r, and C are loaded by the user into the z, y, and x regis- 
ters, respectively. After loading, the user must press "Go to (0)(0)" and "Continue. " 
The program contains a "Pause" in the iteration loop to display p i? |p i+1 - Pjl , and 
C - C^ in the x, y, and z registers, respectively. Keeping the "Pause" button de- 
pressed results in a "Stop" at this display. The final display is p { , I Pf ” Pf- 1 1 > 

C - Cj i in the x, y, and z registers, respectively. Pressing "Continue" after the 
final display sends the calculator through another iteration for added accuracy in p. 
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The following sample cases illustrate the use of the indirect program. The times 
given are the approximate delays between pressing "Continue" and the final display of 
the results. 

Case 4: n = 10 C - C f _ 1 = 0.000 000 562 

r = 7 | Pf - Pf _ 1 1 = 0. 000 000 525 

C = 0.95 p { = 0.849 972 
5 iterations, about 4 seconds 

Case 5: n = 100 C - Cj_j = 0.000 000 019 

r = 75 |p f ~ P f _ x | = 0.000 000 004 

C = 0, 90 p f = 0.797 094 
5 iterations, about 17 seconds 

Case 6: n = 1000 C - C f _ t = 0. 000 005 791 

r = 950 | Pf - p^_ jJ = 0. 000 000 335 

C = 0.95 p f = 0.959 946 
5 iterations, about 32 seconds 


Limitations 

The convergence criterion is based on the increment of p (eq. (9)). Since com- 
parable accuracy does not always exist in C, the error C - Cj_^ is displayed. More 
accuracy in p and, hence, in C can be obtained by pressing "Continue" and going 
through another iteration. 


Error Lights 

The following conditions result in an error light: 

(1) An intermediate p^ + ^ outside the allowable range of 0 to 1. There are two pos- 
sible causes of this condition: (a) C. is incorrectly set to zero when pF is less than 

-99 ^ * 

10 (see eq. (6)); this erroneous Cj is used in equation (5). (b) The derivative is in- 
correct when intermediate steps in computing equation (7) exceed the capacity of the cal- 
culator; this erroneous derivative is used in equation (5). 

Figure 4 shows the allowable pairs of n and r that yield correct solutions for 
C s 0.05. Use of C values less than 0.05 further contracts the allowable region, where- 
as values of C greater than 0. 05 slightly increase the size of the allowable region. For 
n < 300, r can take any value from 2 to n - 1. For 300 < n < 550, there are few re- 
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so 


p = i - (i - C) 1/n (11) 

When r = n, 

C = p n (12) 

SO 

p = C 1 / 11 (13) 

APPLICATIONS 

Several applications of the cumulative binomial function are discussed in the refer- 
ence. In so-called r-out-of-n redundancy problems, the indirect program is useful in de- 
termining the required element reliability p. For binomial sampling plans, either pro- 
gram can be used to determine points on operating characteristic curves. The indirect 
program is useful in determining binomial confidence intervals, whereas the direct pro- 
gram is more conveniently used in hypothesis testing. 


CONCLUDING REMARKS 

The computer programs presented herein for the cumulative binomial function are 
useful for quickly determining the cumulative binomial sum C or, alternatively, the bi- 
nomial parameter p. The probability of exactly r successes among n elements can be 
obtained from the direct program by subtracting the two appropriate cumulative sums; 
that is, the individual r binomial term is obtained by subtraction of the cumulative 
sum over the limits r + 1 to n from the cumulative sum over the limits r to n. 
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